Microwave realization of quasi-one-dimensional systems with correlated disorder 
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A microwave setup for mode-resolved transport measurement in quasi-one-dimensional (quasi- 
1D) structures is presented. We will demonstrate a technique for direct measurement of the Green's 
function of the system. With its help we will investigate quasi- ID structures with various types of 
disorder. We will focus on stratified structures, i.e., structures that are homogeneous perpendicular 
to the direction of wave propagation. In this case the interaction between different channels is 
absent, so wave propagation occurs individually in each open channel. We will apply analytical 
results developed in the theory of one-dimensional (ID) disordered models in order to explain main 
features of the quasi-lD transport. The main focus will be selective transport due to long-range 
correlations in the disorder. In our setup, we can intentionally introduce correlations by changing 
the positions of periodically spaced brass bars of finite thickness. Because of the equivalence of the 
stationary Schrodinger equation and the Helmholtz equation, the result can be directly applied to 
selective electron transport in nanowires, nanostripes, and superlattices. 

PACS numbers: 72.15.Rn, 42.25.Bs, 42.70.Qs 
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I. INTRODUCTION 



Because of progress in nano- and materials science, 
the investigation of (electron) wave propagation in pe- 
riodic one-dimensional (ID) systems has recently at- 
tracted growing attention (see, e.g., Ref. [l] and references 
therein). Systems of particular interest are metamateri- 
als or stacked layer structures, as well as superlattices. 
One of the important problems that is still far from be- 
ing completely understood is the role of disorder due 
to fluctuations in the width of layers or due to varia- 
tions in the material parameters, such as the dielectric 
constant, magnetic permeability, or barrier height (for 
electrons This disorder was seen mainly as an un- 

wanted side effect in experimental realization^. But with 
enhanced ability to control these features, disorder itself 
becomes an interesting candidate for targeted manipula- 
tion of transport properties. In particular, the correla- 
tions in the disorder may exhibit unusual features. One 
of the earliest findings was that short range correlations 
can inhibit localization in low dimensional system a 15 i 16 . 

Later, it was shown, both theoretically-!^— and 
experimentally^ 7 ~— , that specific long-range correla- 
tions can even be used to significantly enhance or 
suppress the wave localization in any desired win- 
dows of frequency. Recent analytical studies extended 
the theoretical predictions from purely ID models to 
quasi-one-dimensional (quasi-lD) randomly stratified 
structure a 23 i 30 , as well as to random, periodic-on-average 
stacked systems with bilayer unit cells^ir— . These stud- 
ies provide the application of correlated disorder to pho- 
tonic crystals, metamatcrials and layered structures, e. g., 



semiconductor superlattices. 

In this paper we explore the validity and applicabil- 
ity of the theoretical results by studying the microwave 
transmission through a quasi-lD waveguide open on both 
ends. It was shown in Ref. that the transport prop- 
erties can be effectively studied by measurement of the 
Green's function, up to a factor describing the antenna 
properties. This method can be considered as an es- 
sential tool to analyze the experimental data. It al- 
lows us to test experimentally the analytical predictions 
from Ref. HH in open waveguides with and without dis- 
order. Furthermore, we extend the experimental confir- 
mation of selective transport to quasi- ID systems that 
emerge because of long-range correlations. Many studies 
on ID systems with correlated disorder involve hundreds 
of scatterer o 16 i 26 i 27 . In present experiments only 26 scat- 
terers are used. We demonstrate that this small number 
is already sufficient to reveal the predicted effect of cor- 
related disorder. 



II. EXPERIMENTAL SETUP AND BASIC 
EQUATIONS 

In our study we use an experimental setup that allows 
us to measure microwave transport through a multimodc 
metallic- wall waveguide with an internal structure of dif- 
ferent kinds (see Fig. [1} . The input (transmitting) and 
output (receiving) antennas are plugged into the waveg- 
uide via slides that can be shifted stepwise by two motors. 
The absorbers are placed on both ends of the waveguide, 
in order to reduce the reflection from the ends, i.e., to 
consider the setup as a finite waveguide with open ends. 
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FIG. 1. (Color online) Not-to-scale sketch of experimental 
setup with top plate removed. Brass-bar inlays can be freely 
placed inside the waveguide. Number of bars varies accord- 
ing to spacing d a . Antennas shown here with parts of the top 
plate can be moved separately along the z direction. Distance 
along the x axis between antennas is L — 1.50 m; total waveg- 
uide length is 2.38 m. Absorbers at the ends of the waveguide 
are omitted. 



The values of the height h, width w and the distance be- 
tween the antennas - the effective length L of the waveg- 
uide - meet the conditions 



h < w <C L 



(2.1) 



that are important for a proper theoretical analysis. 

In an empty rectangular waveguide the electromag- 
netic field can be separated into two components of dif- 
ferent polarization: the TM polarization with transverse 
magnetic field and TE polarization with transverse elec- 
tric fielci^. In our experiments all measurements are 
performed below the cut-off frequency v cu t = c/2h w 
18.75 GHz. In this case only the lowest TE-component, 
having the electric field stretched along the y axis and the 
magnetic field being in the (xz) plane (see the coordinate 
system in Fig. [I}, can propagate. Therefore, our experi- 
mental setup within the frequency interval v < v cnt and 
with the conditions in Eq. (|2.1[) can be considered as a 
quasi-lD open-strip structure, occupying the region 



x < L, 



O^z^w. 



(2.2) 



In order to develop a proper theory for our model, we 
treat the input antenna 1 in (x, z) plane as a point source, 
while the output antenna 2 is an observation point. As 
is known, the electric field of a point source is deter- 
mined by the retarded Green's function Q(\x — x'\;z, z'). 
In Ref. [35| it was shown that for the correct analysis of 
transport one has to take into account the influence of 
coupling between the antennas and the waveguide. As a 
result, the Green's function for the empty waveguide in 



the normal-mode representation reads 



Q{\x - x'\;z, z') = ^ A n sm sin J 



exp (ik n \x — x'\) 
ik n w 



(2-3) 



Here we introduced the factor A n which effectively de- 
scribes the coupling between the antennas and normal 
modes labeled by index n. In our experiment this factor 
A n is regarded as a fitting parameter. Without antenna 
coupling (A n = 1) Eq. 12.31 coincides with the Green's 
function of the two-dimensional Hclmholtz equation com- 
plemented by the Dirichlct boundary conditions at metal- 
lic walls of the waveguide. 

The quantities nn/w and k n are, respectively, the dis- 
crete values of transverse, k z , and longitudinal, k x , wave 
numbers: 

2tt 



k n = \Jk 2 - (im/w) 2 



c 



(nv c 



n= 1,2,3,..., N w . (2.4) 

Here k = 2t:v/c is the total wave number for the elec- 
tromagnetic wave of frequency v. For the modeling 
of electron quasi- ID structures (such as nanowires and 
nanostripes) the wave number k should be regarded as 
the Fermi wave number within the isotropic Fermi-liquid 
model. 

The total number N w of propagating waveguide modes 
(for which k n is real) is determined by the integer part 
[. . .] of the mode parameter kw/n — vjv c , 



N w = {kw/irj = lv/y c \. 



(2.5) 

Correspondingly, the sum in Eq. (|2.3[) runs only over the 
propagating modes with n < N w and ignores the contri- 
bution of evanescent modes for which n > N w . Indeed, 
for evanescent modes the values of k n are purely imagi- 
nary; therefore, they do not contribute to the transport. 
The critical frequency v c , 



c/2w, 



(2.6) 



below which there are only evanescent modes (kw/ir < 1, 
i.e., N w — 0) is called the cut-off frequency. Evidently, 
there is no transport for v < v c . Note that here v c is 
the cut-off frequency for the first mode with n = 1. The 
cut-off frequency for each higher mode with index n 



is defined as u, 



(n) 



III. SCATTERING MATRIX 



It is suitable to define the scattering matrix S nm for 



the whole waveguide 
Fourier transform, 



L) as the twofold sine- 



W Jo \ w J 



w 



Q{L;z,z). 

(3-1) 
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FIG. 2. (Color online) Transmission spectrum \G(L; z, z') 2 
through the empty waveguide for fixed antenna position 
z=z'— 55mm. Vertical lines show the cut-off frequencies ui 
for different modes. 

One can see that in accordance with Eq. (|2.3[) the scatter- 
ing matrix for the empty quasi-lD waveguide [Eq. (|2.2|) ] 
is diagonal in the mode representation. In this case the 
S matrix is proportional to the Green's function for free 
ID propagation with the wave number k n : 

S nm = A n eXp { ^ nL) Q(N w - n)S nm . (3.2) 

Here Q{x) stands for the Heaviside unit-step function, 
where Q(x < 0) = and Q(x ^ 0) = 1. The result 
[Eq. (|3.2[) ] is quite natural, since in an empty waveg- 
uide there are no transitions between the normal modes. 
Therefore, the transport through any channel can be con- 
sidered independently of other channels. 

Below we call the elements S nm of the scattering ma- 
trix [Eq. (|3.Ij) ] the scattering amplitudes. Then, the 
scattering probability T nm from the mode n to mode m, 
the mode-transport coefficient T n (for the scattering of a 
given nth "incoming" mode into all "outgoing" modes) 
and the total transport coefficient T are naturally defined 
as follows: 

N w N w 

^nm — S nm , T n ^ * T n7n ^ T ^ ^ T n . (3.3) 

m— 1 n— 1 

One should stress that although the above definitions 
of T nm , T„ and T correspond to the standard ones known 
in the theory of transport through quasi-ID disordered 
structures (e. g. see Ref . [37h . here the physical meaning of 
all these quantities is different. Indeed, the expression for 
the scattering matrix S n m is originated from the single- 
particle Green's function, not from the Kubo formula, 
which is based on the two-particle Green's function. The 
price for this simplification can be directly seen from 
Eq. (|3.2[) . which has a divergence at values of k n = 0. 
Therefore, the transmission coefficients [Eq. (|3.3[) ] can 
not be normalized. On the other hand, there is no di- 



vergence in the standard expressions for scattering ma- 
trices. However, our point is that this somewhat new 
definition of the scattering matrix, based on the mode 
representation of the Green's function, has an advantage 
when we consider setups similar to the one in our ex- 
periment. In spite of the absence of rigorous expressions 
for the localization length developed in connection with 
this type of scattering matrix, in the following we show 
that the knowledge of theoretical results obtained for ID 
disordered models very much helps us to understand the 
transport properties of the experimental setup. 

For the experimental data analysis, the continuous 
Fourier transform [Eq. (|3.f [) ] is replaced by its discrete 
counterpart, allowing us to compute S nm from Q(L\ z, z') 
measured for different positions z and z' of input and 
output antennas. A typical measurement series with the 
starting points z = z' = 5 mm and step size 10 mm yields 
1 2 individual measurements. An example of the trans- 
mission spectrum obtained in one measurement is shown 
in Fig. [5] It should be kept in mind that in our experi- 
ment the total transport coefficient T is always less than 
I, since the antennas are not perfectly coupled to the 
waveguide. Furthermore, at both ends about half the 
energy escapes toward the open ends. The amount of en- 
ergy lost is not accessible in the present setup. Moreover, 
the absorption inside the waveguide reduces the value of 
T further. Practically, in our experiment the total trans- 
port coefficient T never exceeds the value T = 0.17. 

The Fourier transform [Eq. (|3.ip ] can be performed for 
every single frequency v in the spectra yielding the scat- 
tering matrix S nm as a function of v. An example of the 
scattering probability T nm = IS^ml 2 for an empty waveg- 
uide and fixed value v — f f GHz is shown in Fig. [3] At 
this frequency there are seven propagating modes. Since 
there are no scatterers in the waveguide, there should 
be no transitions between different propagating modes, 
therefore, the matrix T nm is expected to be diagonal in 
accordance with Eq. (|3.2p . Experimentally, the total con- 
tribution of the off-diagonal elements in comparison with 
that of diagonal elements is found to be always less than 
5 % (typically much smaller) . The main source of these 
off-diagonal terms is the perturbation caused by anten- 
nas. We always observe that the highest propagating 
mode gives the main contribution to the transport, as in 
Fig.[3J where the term 15771 2 is the largest one. This fact 
is in complete correspondence with the expectation from 
the theoretical predictions in the next section. 



IV. EMPTY WAVEGUIDE 

Let us analyze the transport through various modes as 
a function of the wave frequency. As an example, the 
measured frequency dependence of the mode-transport 
coefficient [Eq. (|3.3|) ] for the third T 3 (i/) and the fourth 
T±(v) modes is depicted in Fig. @] by solid curves. There 
are several interesting features in this figure. First of 
all, there is a pronounced spike in T%{v) at the cut-off 
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FIG. 3. (Color online) Modulus square of the scattering ma- 
trix elements T nm = \S„m\ 2 at v = 11 GHz with seven prop- 
agating modes. 

(3) 

frequency v c = 4.5 GHz, at which the third mode opens 

(3) 

and starts to propagate. The experimental value of v c 

(3) 

coincides very well with the theoretical prediction, v\ = 
3c/2w [see Eq. (|2.6p and comments below it]. Above this 
critical value the coefficient T3 decreases for increasing 
frequency. 

The frequency dependence of T^(y) looks quite similar 
to that of T$(y), apart from the shift of the spectrum 
by the value v c = c/2w = 1.5 GHz. This is because the 
propagation of the fourth mode emerges above the cut- 
off frequency v[ 4 ^ = 4c/2w = 6.0 GHz. In general, all 
mode-transport coefficients T n (v) have somewhat simi- 
lar frequency dependence. Namely, they manifest the 
corresponding spikes at the corresponding cut-off values, 

i>c = nc/2w, and decrease with a further frequency in- 
crease. 

In order to understand such specific frequency depen- 
dence of the coefficients T n (u), one has to use Eqs. p.2[) 
and (|2.4[) . Since in an empty waveguide all propagating 
modes are independent of each other, the mode-transport 
coefficient reads 

T n {v) = \S nn \ 2 = ^r^6(N w - n) 

r 2 \A I 2 

= 1fi 2 ( 2 2 2 M "-nv c ). (4.1) 

One can easily achieve a nice correspondence of this the- 
oretical result to the experimental data by a simple fit 
of the factor A n . To do this, we must also take into ac- 
count an evident assumption that the wave number fc, 
(i.e., the frequency v) has a small imaginary part Irak 



(or Im v) , which emerges because of an absorption in the 
waveguide. 

Typically, the absorption is frequency dependent, but 
in the analysis it is sufficient to treat it as a constant 
for each mode, i.e., to regard Imfc as a function of the 
mode index n only. In Fig. 0] the theoretical expression 
(|4.ip is displayed by dotted curves. Note, however, that 
the fast oscillations that are clearly seen in the experi- 
mental curves, are not described by our theory. These 
oscillations are caused by multiple wave reflections from 
the waveguide ends (occurring because of non-perfect ab- 
sorbers) and from the antennas. 

The frequency dependence experimentally observed for 
the coefficient T n {v) is explained by the presence of the 
squared modulus of longitudinal wave number k n in the 
denominator of Eq. (|4.1[) . Apart from the frequency de- 
pendence, this term, \k n \ 2 = ^-(v 2 —n 2 v 2 ) 1 provides the 
global dependence of T n {v) on the mode index n. One 
can see that at fixed frequency v the higher the mode 
index n, the larger the coefficient T n (y). That is why the 
highest propagating mode with n = N w always gives the 
main contribution to the total transport coefficient T, es- 
pecially, in the vicinity of its cut-off, v^ Nw ' = N w c/2w. 
This effect is displayed in Fig. [31 



V. PERIODIC STRUCTURE 

The features of the empty waveguide discussed above, 
were explored as a necessary precondition for our fur- 
ther analysis. As a second step, we studied a periodic 
structure without any disorder. To this end, the brass 
bars of thickness df, = 8 mm measured along the x axis, 
were inserted into the waveguide as is shown in Fig. [T] 
The constant spacing d a between neighboring bars was 
chosen as d a = 17 mm and d a = 32 mm for two differ- 
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cnt realizations, yielding a period d = d a + d = 25 mm 
for the first realization with 42 bars, and d = 40 mm for 
the second one with 26 bars. The measurements were 
performed as before, giving rise to the scattering matrix 
S nm of the periodic setup. Both periodic structures gave 
similar results; therefore, in what follows, we restrict the 
discussion to the second case only. It would be inter- 
esting to investigate the behavior of the bands and the 
correlation gaps with decreasing number of scattcrers N, 
but this is beyond the scope of the present paper. 

Since the periodic structure is arranged along the wave 
propagation (along the x axis) and the bars are homo- 
geneous along the z axis, the propagating modes remain 
independent as in the empty waveguide. Indeed, the scat- 
tering matrix [Eq. (|3.1|) ] was experimentally found to be 
diagonal within the same accuracy as before. As a con- 
sequence, the spectrum of the mode-transport coefficient 
T n (v) = \S nn \ 2 looks similar for all propagating modes, 
provided the shift of spectrum by particular cut-off fre- 
quencies, = nc/2w is done. This means that T n (v) is 
convenient to represent as a function of the correspond- 
ing longitudinal wave number k x = k n [see Eq. (|2.4[) ]. 
Fig. EJa) shows the semilogarithmic experimental plot of 
the mode-transport coefficient T3QC3) for the third prop- 
agating mode. One can clearly see a band structure in 
the dependence T^fca), as is expected for a periodic ar- 
rangement. Note that the experimental transmission in 
the gaps is smaller than in the bands over several orders 
of magnitude. 

From the physical point of view our setup with peri- 
odic filling of brass bars represents a set of N w indepen- 
dent ID periodic arrays of two alternating a and b layers. 
Note, however, that each unit (a, b) cell of every nth ar- 
ray has a quite complicated structure. While the first 
a-laycr is a trivially homogeneous air spacing of thick- 
ness d a with the wave number k n , the second &-layer of 
thickness <4 is not homogeneous along the y-axis. Indeed, 
since the height (5 mm) of the brass bar is smaller than 
the waveguide height (h = 8 mm), there is an air slot be- 
tween the metallic bar and the waveguide top plate, see 
Fig. [1] For this reason, in our analysis we shall charac- 
terize the 6-layers by the wave number k = n k ni thus, 
introducing an effective refractive index n& that should 
be fixed by the experimental data. In this way, in or- 
der to describe the transmission band structure for every 
nth array of bilayers (for every nth propagating mode) 
we employ the effective Kronig-Pcnncy model with the 
common dispersion relation for the Bloch wave number 
K (see, e.g., Refs. [H andl38h. 

cos(nd) = cos(k n d a ) cos(n k n db) 

~a + sin(k n d a )sm(n b k n d b ). (5.1) 

For the standard bilayer models the factor a+ is the 
arithmetic average of the impedance ratio and its inverse 
value, a + = (Z a /Z + Z /Z a )/2. Since the impedance Z 
of Mayers is unknown for our system, the factor a + is 
regarded as the second fitting parameter. Note that a+ 
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FIG. 5. (Color online) Band structure for (a) periodic (b) 
white noise and (c) correlated disorder arrangement with 
{da(j)) — d a = 32mm. Shaded regions indicate the bands 
restricted by phenomenologically obtained band edges (ver- 
tical lines). Dashed vertical line indicates the Fabry-Perot 
resonance (see Section [VI All . 

equals one for the case of two basic layers a and b which 
are optically matched (Z a = Z ). 

The solution of Eq. (|5 . 1 [) gives the Bloch wave number 
re = n(k n ) as a function of the mode wave number k n 
(or, as a function of the wave frequency v and the mode 
index n) . The real values of re are achieved for such k n for 
which I cos(recZ) < 1, thus, resulting in the emergence of 
spectral transmission bands. Otherwise, for k n for which 
|cos(red)| > 1, the Bloch wave number re turns out to 
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FIG. 6. (Color online) The r.h.s. of the dispersion relation 
(|5.ip (solid line) for the best fit of the parameters. The spac- 
ing is d a = 32 mm with fitting parameters a+ = 1.6 and 
rib = 1.42. Full circles correspond to the band edges extracted 
from Fig. Oa). 

be imaginary leading to the emergence of spectral gaps. 
The dispersion relation (|5.1[) gives the band edges at the 
points where cos{nd) = ±1. 

As a first step, we manually extracted the positions 
of the band edges from the measured spectrum of the 
coefficient T n {y) = \S nn \f , see vertical lines in Fig. [5] for 
these positions. After, we fitted the dispersion relation 
(|5.1[) by changing the values of a+ and rib, and keeping 
the crossing points of its r.h.s. with the horizontal lines 
coa{nd) = ±1, in the correspondence with the found band 
edges. The result is displayed in Fig. [B] for a + = 1.6 and 
n b = 1.42. The dark dots in the figure show the extracted 
band edges. 

VI. POSITIONAL DISORDER 

In order to study the effect of disorder, we shift each 
brass bar by a random value along the waveguide length 
(the x-axis), keeping the average distance between the 
bars constant. Since the disorder is a function of the 
x-coordinate only, the propagating modes remain inde- 
pendent as in the considered above arrangement s 23 ' 30 . In 
such a manner, we can consider every nth channel as an 
independent one-dimensional array of two alternating a 
(air spacing) and b (brass bar) layers. Theoretically, the 
positional disorder is incorporated via the random thick- 
ness of a-layers only, 

d a (j) = d a + Qa(j) , (d a (j)) = d a , (6.1) 

while the thickness d of b layers is constant. The integer 
j = 1, 2, 3, . . . , Nd, with Nd = \L/d\ enumerates the unit 
(a, b) cells, d a is the average thickness of the a layer; and 
Qa(J) stands for small random variations. Note that in 
the absence of disorder (g a = 0) the set of ID bilayer 
arrays represents the periodic structure with the period 
d = d a + d that was considered in previous section. 



The random sequence g a (j) is statistically homoge- 
neous with the zero average and given variance <r 2 , 

(&.(?')> =0, (lfiV)) = o*. (6.2) 

The binary correlator is defined by 

(Qa(j)Qa(j')} = crlK a (j-j'). (6.3) 

Averaging (...) is performed over the whole array of lay- 
ers or over different ensembles, which are assumed to 
be equivalent. As one can see, the two-point correlator 
K a (j — j') is normalized to unity, K a (0) = 1. 

If the positional disorder is weak (fc 2 cr 2 <C 1), all trans- 
port properties are determined entirely by the random- 
ness power spectrum K, a {k), 

oo 

K a {k) = 1 + 2 K a (r) cos(fcr), (6.4a) 

r—l 

i r 

K a(r) = - dkK, a {k)cos{kr). (6.4b) 
^ Jo 

Since the correlator K a (r) is a real and even func- 
tion of the difference r = j — j' between cell indices, 
its Fourier transform lC a (k) is a real and even function 
of the dimcnsionless wave number k. According to the 
Wicncr-Khinchin theorem, the power spectrum K. a {k) is 
non-negative for any real sequence g a (j)- 

The transport through any ID disordered structure is 
governed by the phenomenon of Anderson localization 
(see, e.g., Ref . l39h . Its principal concept is that all trans- 
port characteristics depend only on the ratio between the 
system length L and localization length L\ oc . Thus, the 
localization length is the key quantity that controls the 
transport in a ID geometry The analytical expression for 
the localization length in the periodic-on-average model 
with bilayer structure was derived in Ref. [U Specifi- 
cally, the case of weak disorder in the thickness of only 
one kind of layers was analyzed, which is close to our 
experimental setup. 

However, our model is essentially quasi-one- 
dimensional. On the other hand, in Refs. [23| and 
l30l it was shown that in quasi- ID waveguides with 
random stratification along the propagation, each of 
the (independent) modes can be associated with their 
own length Li oc (k n ). The latter can be obtained by 
a direct substitution of specific mode characteristics 
into the general expression for the localization length 
of ID random structures. Applying the results from 
Refs. HH, US and HH, one can readily get the following 
expression for the inverse localization length, associated 
with the nth propagating mode of our setup, 

L iol( k n) = r, , ■ 2/ ° £ a (2KcQsin 2 (n b fc„(4). (6.5) 
2d sin (Kd) 

Here a 2 _ is the so-called mismatching factor. For stan- 
dard bilayer models a_ is related to a + in the disper- 
sion relation [Eq. (|5.1|) ] by the equality a\ — o?_ = 1. 
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FIG. 7. Ratio L/L [oc (k n ) formed by Eq. (|6.5j) for white- 
noise disorder (solid curve) and for correlated disorder (dot- 
ted curve). Parameters are d a = 32mm, a a = 2.21 mm, a 2 _ = 
1.56. Horizontal line indicates the equality L = Li oc (k n ). 



Therefore, for the perfect matching for which ct + = 1, we 
have a_ = resulting in a divergence of the localization 
length, Lioc = oo. Since for our system a+ = 1.6, the 
mismatching factor is a 2 _ = 1.56. Note that Eq. (|6.5p is 
a particular case of a more general expression derived in 
Ref. |32| for the case of both a and b positional disorders 
(when thicknesses of both a and b layers are randomly 
perturbed). 



A. White-noise disorder 

In Fig. [7] the inverse localization length L^ c (k n ) of 
the nth propagating mode, measured in the units of the 
waveguide length L = 1.5 m, is plotted against k n . The 
26 displacements g a (j) were drawn randomly from the in- 
terval -5 mm to 5 mm, while the variance was constrained 
to <r 2 = 4.9 mm 2 . The ratio of a a and the spacing d a is 
about 7%. 

Now we start with the discussion of experimental 
data for the case of uncorrelated (white-noise) disor- 
der. In this case the binary correlator is K a (j — j') — 
djji , and consequently, the power spectrum is constant, 
K-a{k) = 1. The corresponding ratio L / Li oc {k n ) is pre- 
sented in Fig. [7] by solid curve. One can see that within 
the first spectral band the inverse localization length 
monotonously increases from zero to infinity dependen- 
ing on k n . Within any other Bloch band the mode- 
localization length Li oc {k n ) vanishes at the band edges 
(where Hid — 0,7r) because of the term sin 2 (re<f) in the 
denominator of Eq. (|6.5|) . 

It is important that the mode-localization length 
Lio C (k n ) exhibits the Fabry-Perot resonances associated 
with multiple wave-reflections inside every b layer from 
its interfaces. As is known, the resonances appear when 
the constant width db is equal to an integer multiple of 



half of the wavelength inside the layer: 

k n — Trs/ribdb , s = 1, 2, 3, . . . . (6.6) 

At the resonances the factor sin 2 (n£,fc„<4) in Eq. (|6.5[) 
vanishes, resulting in resonance divergence of the local- 
ization length and, consequently, in suppression of the 
localization (the nth mode becomes fully transparent). 
Remarkably, the Fabry-Perot resonance should be quite 
broad, since it is caused by the vanishing of a smooth 
trigonometric function. This effect was analyzed both 
theoretically and experimentally in Refs.[33| andHJ. Since 
in our experimental setup rib = 1.42 and db = 8 mm, 
the first (s = 1) Fabry-Perot resonance is expected at 
k n = 280.5 m _1 and this is what is observed in the fourth 
transmission band in Fig. [7] 

In accordance with single-parameter scaling, there are 
only two transport regimes in the ID Anderson localiza- 
tion: the ballistic (L <C Li oc ) and localized (Li oc <C L) 
regimes. In the ballistic regime a sample is practically 
transparent, and in the localized transport a ID disor- 
dered structure almost perfectly reflects classical or quan- 
tum waves. 

In Fig. [7] the horizontal straight line L/Li oc (k n ) = 1 
separates the upper region, L/Li oc (k n ) > 1 (where the 
nth propagating mode is localized and therefore closed), 
from the lower one, L / Li oc {k n ) < 1 (where the nth propa- 
gating mode is transparent). As one can see from Fig. [71 
in our experiments the first and fourth spectral bands 
should be almost completely transparent. Note that the 
transparency of the first band is due to a very weak in- 
fluence of the disorder, while the transparency of fourth 
band is caused by the broad Fabry-Perot resonance. 

The predicted features of the system are clearly dis- 
played in Figs. [5] for the experimentally measured mode- 
transport coefficient T^k^). The comparison of Figs. [5] 
(a) and (b) shows that the introduced disorder does not 
reduce transmission in the first band, where k n < 50 m _1 . 
In the fourth band between 240 m _1 and 290 m -1 again 
there is only a weak influence of disorder. Only in the 
second, third and fifth bands the compositional disorder 
reduces the mode transport by 2-3 orders of magnitude. 



B. Correlated disorder 

As one can see from Eq. (|6.5p . the inverse localization 
length LT~z(k n ) follows the profile of the power spectrum 
]C a (2nd). This fact is of special interest in view of exper- 
imental realizations of positional disorder with specific 
long-range correlations. In particular, one can artificially 
construct an array of random bilayers with such a power 
spectrum that vanishes abruptly within prescribed inter- 
vals of the mode wave number k n , resulting in the diver- 
gence of the mode-localization length. Thus, by a proper 
design of the positional disorder, one can create quasi- 
1D bilayer structure with a selective mode transmission 
T n (kn) in dependence on the frequency v. 
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FIG. 8. Power spectrum /C a (2Kc0 versus dimensionless Bloch 
phase nd/ir for correlated (dotted curve) and uncorrelated 
(solid curve) random sequences. The factor K. a (2nd) is calcu- 
lated from Eq. (|6.4[) taking only the first 13 terms in the sum 
that corresponds to an arrangement with Nj = = 26 

unit (a, b) cells with (d a (j)) — d a — 32 mm. 



To do this, we perturbed all periodically spaced = 
26 brass bars by random values Q a (j) (j = 1,2,3,..., N d ) 
correlated in such a way that the power spectrum 
K, a {2nd) is zero everywhere, except for narrow regions 
around nd/n = 0.3 and 0.7. 

We can solve the problem of generating the random se- 
quences with prescribed short- or long-range correlators 
by employing a widely used convolution method origi- 
nally proposed by Rice^S. The details of various appli- 
cations of this method can be found in Refs. Ii~7l - fl9l I2TI — 
HI, [M H3, US EHH With the use of this method, the 
resulting curve for K, a (2nd) has the form shown in Fig. [5] 
For comparison, the power spectrum for an uncorrelated 
random sequence is also shown. The inverse dimension- 
less localization length L/Li oc (k n ) that corresponds to 
the specified correlated disorder, is depicted in Fig. [7] by 
the dotted curve. In accordance with this curve, one 
should expect the appearance of correlation gaps, i.e., 
additional gaps due to correlated disorder (namely, one 
correlation gap in the first band and two in the second 
and third transmission bands). In the first band the first 
possible correlation gap at k„ ~ 25 m" 1 remains in the 
ballistic regime and can not be observed. The same hap- 
pens in the fourth band, where only the first correla- 
tion gap can be observed. Because of the overlap with 
the Fabry-Perot resonance the second correlation gap is 
shifted to the ballistic regime. 

Fig. [He) shows the mode-transmission through the 
correlated random arrangement. Globally, the spectral 
structure looks similar to that for the uncorrelated ran- 
dom setup. Indeed, we can see that the transport de- 
creases with increase of k n , with an exception of the 
fourth band containing the Fabry-Perot resonance. In 
the first band the mode-transport coefficient T n (k n ) is 
comparable to that for both regular and uncorrelated 
random setup: Figs.[SIa) and (b). This is because in this 
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FIG. 9. (Color online) Comparison of the mode-transport 
through uncorrelated (solid curve) and correlated (dotted 
curve) random arrangements. Upper figure shows the second 
band; lower figure presents the fourth (Fabry-Perot) band. 
Vertical dotted line indicates expected correlation gaps for 
Hid/ 7t = 0.3, 0.7. Vertical dashed line indicates the position of 
the Fabry-Perot resonance at k n = 280.5 m" 1 . 



band the localization length is larger than the system 
size, L <C Li oc (k n ). Only in the region of the correlation 
gap at k n = 40 m" 1 does the localization length shrink, 
such that Li oc (k n ) L, resulting in an very weak, but 
recognizable correlation gap. This correlation gap is an 
interesting special case, because it is in the regime inter- 
mediate between ballistic and localized. 

The second band exhibits two correlation gaps that 
are clearly seen, in spite of the fact that in this band the 
transport is reduced by two orders of magnitude. In the 
third spectral band, at about k n = 200 m" 1 , the trans- 
port already decreased to noise level. However, in the 
fourth band containing the Fabry-Perot resonance, the 
transport is quite large and comparable to that for the 
regular arrangement - apart from one correlation gap. 
Finally, within the fifth band the mode-transport coef- 
ficient is reduced by randomness so strongly that it is 
not possible to reveal any difference between correlated 
and uncorrelated disorder. All these properties are in a 
complete agreement with the predictions provided by our 
theoretical model. 

In order to compare the effects of uncorrelated and cor- 
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related disorder in more detail, we present Fig. [9j which 
illustrates the second and fourth spectral bands of the 
mode-transport coefficient T n (k n ) for the third propa- 
gating mode. The expected positions of the correlation 
gaps at nd/ir = 0.3, 0.7 are indicated by vertical dotted 
lines, and the position of the Fabry-Perot resonance is 
given by the vertical dashed line. 

In the second spectral band (see upper part of Fig. [9}, 
the correlation gaps are easily recognized, showing the 
variation of the amplitude over five orders of magnitude. 
There is good agreement between the theoretically pre- 
dicted and experimentally found positions of these cor- 
relation gaps. Between the correlation gaps, i.e., in the 
vicinity of k x = 110 m , the mode-transport coefficient 
for the correlated disorder is two orders of magnitude 
higher than for uncorrelated disorder. Again, this is 
not surprising, since the power spectrum /C a (fc) ~ 1 for 
uncorrelated disorder, whereas in the case of correlated 
disorder it is close to zero outside the correlation gaps. 
Remember that the mode-localization length [Eq. (|6.5p ] 
tends to infinity as JC a {k) — > 0. 

In the fourth spectral band containing the Fabry-Perot 
resonance (lower part of Fig. |9]) only the first correla- 
tion gap is found. The second correlation gap is at a 
short distance from the Fabry-Perot resonance and there- 
fore, is suppressed. Indeed, here because of the reso- 
nance the mode-localization length is still quite large, 
L <C Li oc (k n ), in order to contribute to transmission 
of the waveguide of finite length, see also Fig. [7] Only 
for the first correlation gap does the mode-localization 
length turn out to be sufficiently small, Li oc (k n ) <C L, 
to observe this gap. Note that for white-noise disorder 
the fourth band is almost fully transparent exclusively 
because of the Fabry-Perot resonance. Nonetheless, the 
specific correlations in the positional disorder of the brass 
bars can induce a correlation gap provided that the dis- 
tance between the resonance and the corresponding max- 
imum of the power spectrum IC a (2Kd), are sufficiently 
long. 

Finally, we would like to stress that all peculiar effects 
discussed above, are experimentally observed for a quite 
short quasi-lD system, with as few as Nd = 26 brass 
bars. This fact confirms the expectation that the analyt- 
ical expression for the localization length, formally ob- 
tained for infinite systems, is actually applicable for real 
experimental situations with a relatively small number of 
scatterers. 



VII. CONCLUSIONS 

We studied transport through a quasi-lD microwave 
waveguide with periodic and randomized scattering po- 
tentials. The potentials were stratified, i.e., homogeneous 
along the width of the waveguide. In this case the prop- 
agation of waves occurs independently along each open 
channel, thus allowing one to apply the theory of ID 
transport to any partial channel and make conclusions 



about both partial and total transmission. 

In the experimental setup the incoming and outgoing 
waves are induced by pointlike antennas. Therefore, one 
of the problems was to understand the influence of the 
coupling of antennas to the waveguides, as well as the 
role of reflection from the antennas and the open ends. 
Also, the effect of absorption in the walls influenced the 
transport characteristics. We have tried to incorporate 
all these effects in our phenomenological theory aiming 
to explain the observed data. 

Main attention was paid to selective transport due to 
long-range correlations in the disorder. We generated the 
disorder by intentionally shifting the position of bulk 
scatterers. In order to quantify our experimental setup, 
we started with the case of an open waveguide with- 
out any scatterers. In this case we have confirmed that 
the coupling between various channels due to antennas is 
small and can be neglected when compared to the theo- 
retical predictions. Then we studied scattering through 
the waveguide with periodically inserted scatterers cho- 
sen in the form of finite-thickness brass bars. With such 
an arrangement, the transmission spectrum was found to 
have a band structure, as expected. However, in order 
to analytically describe the structure of the spectrum, 
one must create a phenomenological model that properly 
corresponds to the experimental setup. 

The specific problem we encountered was that the pe- 
riodic brass bars inserted in the waveguide cannot sim- 
ply be considered as finite-thickness barriers with certain 
refractive indices. They were made of metal, and in or- 
der to allow propagation they have a finite height that 
is less than the vertical dimension of the waveguide. In 
this case it impractically to rigorously create a relatively 
simple and adequate theoretical model, and instead we 
decided to use the standard dispersion relation known 
in the theory of ID models, by introducing two fitting 
parameters. In this way we were able to give an ana- 
lytical expression for the positions of transmission bands 
and spectral gaps. This expression helped us describe 
another experiment, in which the positions of brass bars 
were slightly perturbed, both in a random manner and 
with specific long-range correlations. 

For the case of uncorrelated disorder, we obtained ex- 
perimental data that can be described by an analytical 
expression for the localization length. For any kind of 
correlation, the corresponding expression was derived in 
Refs. [H, [30M32I and one of our interests was to check how 
well our experimental data correspond to the predictions 
of the theory. It should be stressed that this problem is 
far from trivial, since the expression for the localization 
length was developed for the infinite number of scatter- 
ers and for a weak disorder. In our case, the number of 
scattering bars Nd = 26 was relatively small, and many 
other experimental imperfections (absorbtion, non per- 
fect coupling, etc.) are not taken into account in the 
expression for the localization length. However, the ex- 
perimental data were found to correspond well with our 
expectations following from the theory. In particular, a 



10 



quite unusual effect of Fabry-Perot resonances was ob- 
served in experimental data, which is fully explained by 
the theory. 

The Fabry-Perot resonances were predicted to play a 
specific role in the transport (see discussion in Refs. [HI 
and HH) . Namely, if one such resonance emerges inside 
the transmission band, the transmission coefficient turns 
out to be very large as compared with the bands with- 
out such a resonance. The nontrivial point is that the 
resonance effect occurs not only for a specific frequency, 
but in a large interval of frequencies inside the band. In 
our experiment we have clearly observed the Fabry-Perot 
resonance in the fourth band. The transport in this band 
would be negligible otherwise. 

Another specific interest was to check whether in our 
setup with a finite number of brass-bar scatterers one 
can observe the effect of long-range correlations embed- 
ded intentionally in the positions of bars. According 
to the theory, with specific long-range correlations one 
can construct such potentials for which the localization 
length (obtained in the first order perturbation theory) 
diverges (or very small) in finite windows of frequency In 
these windows one can observe either strongly enhanced 



or strongly suppressed transport, compared to the white- 
noise disorder case. In order to observe experimentally 
the effect of suppression of the transmission, we created 
disordered perturbations of the positions of brass bars, 
in accordance with theoretical expressions for long-range 
correlations. Remarkably, our experimental data clearly 
manifest the emergence of new gaps inside the transmis- 
sion bands, in good correspondence with analytical ex- 
pressions. In particular, these gaps were found to arise 
at the positions that were predicted theoretically. 
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